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Abstract 

We study some phenomenological models in a matrix model corresponding to 
the IIB matrix model compactified on a six-dimensional torus with magnetic 
fluxes. Extending our previous works, we examine a wider class of models: a 
Pati-Salam-like model with a gauge group U(4) x 11^(2) x \Jji{2), and models 
where the gauge group U(4) is broken down to Uc(3) x U(l) and/or the \Jr{2) 
is broken to U(l)^. We find all the matrix configurations that yield matter 
content of all the phenomenological models whose gauge group is a subgroup 
of U(8). We then estimate semiclassically a probability distribution for the 
appearance of the phenomenological models. 
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1 Introduction 



The standard model (SM) of particle physics agrees well with experiments. 
While phenomenological models beyond the SM will be explored in (near-)future 
experiments, some theoretical guides may be helpful. It may also be important 
to reconsider why the SM is so successful and what we should ask nature next. 

On the other hand, the SM is unsatisfactory as a final theory, and the string 
theory is expected to be an ultimate theory including gravity. Phenomenological 
models inspired by the string theories have been studied extensively (see, for 
instance, ref. [Ij). However, a serious problem in the string theory is that it has 
too many vacua. 

A candidate for a nonperturbative formulation of string theory is matrix 
models (MM) [2l[3lll]. Since the MM has a definite action and measure, we can, 
in principle, dynamically compare the string vacua, and calculate everything, 
such as spacetime dimensions, gauge groups, and matter contents. Indeed, 
spacetime structures have been analyzed intensively, and four-dimensionality 
seems to be preferred in the IIB matrix model [5l|6l[7j. Then, assuming that our 
spacetime is obtained, we will consider how the SM and some phenomenological 
models appear form the MM, and estimate probability distribution of their 
appearance. Such phenomenological studies in the MM may give us a guide for 
exploring phenomenological models beyond the SM. 

An important ingredient of the SM is the chirality of fermions. We usually 
obtain a chiral spectrum on our spacetime by introducing nontrivial topologies, 
which then give chiral zero modes, in the extra dimensions: Euler characteristics 
of compactified manifolds, special boundary conditions at orbifold singularities, 
the intersection numbers of D-branes, etc., give nontrivial topologies. Also from 
the MM, chiral fermions and the SM matter content were obtained by consid- 
ering toroidal compactifications with magnetic fluxes [SJ [S] and intersecting 
D-branes |loj^ . The former is similar to the magnetized D-branes wrapping on 
a torus, which are T-dual to the latter. 

In this paper, we will study phenomenological models in the MM compact- 
ified on a torus, extending our previous works OS]- We here examine a wider 
class of phenomenological models than in ref. [9]. In order to embed the SM 
fermions with three generations into matrices in our formulation, the gauge 
group must be a subgroup of U(8) or larger groups. We then find all the matrix 
configurations that yield all the phenomenological models whose gauge group 
is a subgroup of U(8). Exhausting all the solutions is necessary for studying a 

^ Studies based on fuzzy spheres were given in [111 1121 [13] . MM's for orbifolds and orien- 
tifolds were studied in [141 115[ . 
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semiclassical analysis. 

We here note that even if a matter field is massless at the tree level in the 
MM, which might be interpreted as phenomena at the Planck scale, it would 
obtain a mass through quantum corrections at low energies. While masses of 
the gauge field and fermionic field can be protected by the gauge and the chiral 
symmetry, respectively, a mass of a scalar field is difficult to protect, which is 
well-known as the naturalness or the hierarchy problem. Thus, although the 
matrices provide scalar fields with the same representation under the gauge 
group as the Higgs field, it is difficult to keep them massless. Then, we will 
first find matrix configurations that provide the gauge and the fermionic fields, 
assuming that the electroweak symmetry breaking is caused by those Higgs 
candidates or some other mechanisms. We will next assume situations where the 
Higgs mass is protected by the super symmetry, and find matrix configurations 
that provide the Higgsino fields as well. 

We then study the dynamics of MM semiclassically, and estimate a prob- 
ability distribution for the appearance of the phenomenological model^. The 
fact that we could perform these analyses is an advantage of Mivjfl. On the other 
hand, there remain some important issues in the formulations of MM, such as 
the relations between MM and string theory, interpretations of spacetime and 
matter in the matrices, and how to take large-N limits. We can expect reversely 
that these phenomenological studies may give some hints for those problems in 
MM. 

In section [21 we briefiy review a formulation of topological configurations 
on a torus, while a similar review was given in ref. [9]. We then find matrix 
configurations that provide the phenomenological models without the Higgsino 
fields in section [3l and those with the Higgsino fields in section [H In section 
El we estimate the probability distribution over the phenomenological models. 
Section [6] is devoted to conclusions and discussions. In appendix |Al detailed 
calculations for determining the fluxes are shown. 



2 Review of topological configurations on a torus 

We begin with a brief review of the IIB MM [3]. Its action has a simple form 
5iiBMM = tr (j[Am,An][A^',A^] + hr^'iAM,^]) , (2.1) 

5lIBMM V'* ^ / 

^ Probability distributions over the string landscape were estimated, for instance, by 
counting the number of flux vacua [T7], and by considering cosmological evolutions |18) . 
^ A related work is given in the MM of noncritical strings [19) . 
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where Am and ^ are NxN Hermitian matrices. They are also a ten-dimensional 
vector and a Majorana-Weyl spinor, respectively. Performing a kind of func- 
tional integration 

j clA e^^"BMM (2.2) 

as a statistical system, and taking a suitable large- limit, one can obtain a 
nonperturbative formulation of string theory. Since the measure as well as the 
action is defined definitely, we can calculate everything in principle. 

A notable feature is that both spacetime and matter emerge from the matri- 
ces. Spacetime can be interpreted as an eigenvalue distribution of the bosonic 
matrices Am, and its structures have been analyzed dynamically [3 El [7]. Ar- 
guments for the local fields on it were given [20]. It was also shown that non- 
commutative (NC) space and matter fields on it are described in MM rather 
elegantly [IB 122^ . 

We now assume compactifications to x with carrying nontrivial 
topologies. For concreteness, we consider toroidal compactifications of x 
T^. Toroidal compactifications were studied in Hermitian matrices |24[ [2T] 
and in unitary matrices [25]. We here use a finite-unitary-matrix formulation 
for NC tori. It can be defined by the twisted Eguchi-Kawai model |26i [27] 
(see, for instance, ref. [28]). Note, however, that such details of formulations 
are not relevant for studying phenomenological models. Any compactifications 
with nontrivial topologies can work as well. We then consider background 
configurations corresponding to 

Afj, ~ (g) 1 , 

e*^' ~ l^Vi , (2.3) 

with /i = 0, . . . , 3 and z = 4, . . . , 9. Matrices Vi represent T^, while represent 
our spacetime M^. We do not specify details of our spacetime x^, i.e., whether 
it is commutative or NC, compactified or not. A more precise correspondence 
between the IIB MM and the unitary MM will be given in section [5j 

We now focus on Vi in (j2.3p . i.e., NC with nontrivial topologies. It is 
known that nontrivial topological sectors are defined by the so-called modules 
in NC geometries (see, for instance, ref. [29]). In the MM formulations, the 
modules are defined by imposing twisted boundary conditions on the matrices 
[28\ 130] . In fact, each theory with twisted boundary conditions yields a sin- 
gle topological sector specified by the boundary conditions [311 [32], while in 
ordinary gauge theories on commutative spaces, a theory, for instance, with 

* It was also shown that curved spaces can be described by interpreting the matrices as 
differential operators [23] . 
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periodic boundary conditions, provides all the topological sectors. However, 
since we now want to derive everything from the IIB MM, those topological 
features of NC gauge theories are not desirable. We thus introduce nontriv- 
ial topological sectors by background matrix configurations, not by imposing 
twisted boundary conditions by hand. Nontrivial topologies can be given by 
block-diagonal matrices [8]. We then consider the following configurations: 



3+j 



5+j 



Vn 



7+j 



V 



VphJ 



VphJ 



(2.4) 



with j = 1,2. The number of blocks is denoted by h. Each block is a tensor 



product of four factors. The first three factors each represent T'^ of 



T2 X 



T'^xT'^, and the last factor provides a gauge group structure. The configuration 
(12. 4|) gives the gauge group U(p^) x U(p^) x • • • x U(p''). 

The matrices ■ with a = 1, . . . ,h and / = 1, 2, 3 in ()2.4p are defined by 
using the Morita equivalence. For details, see, for instance, ref. [29 t [28 t \30 \ [8]. 
r"^- are U(n") matrices that satisfy the 't Hooft-Weyl algebra 



— O-wi !- 

6 ' J- L2^ Ll 



(2.5) 



where the integers mf and nf are specified by an integer as 



-si + kqf , 



n, 



Ni - 2riqf , 



(2.6) 



for each a and I. The integers A'^;, r^, si and ki for each / specify the original 
torus of the Morita equivalence for each T'^. Equations ()2.6p can be inverted 



1 = 2rimf + kinf , qf = Nimf + sinf . 



(2.7) 



^ In the notations of p and q in ref. [291 1281 130| . the present case corresponds to p = 1, 
q — qf , and po = p", and thus p — p"^ and g = p"?;". Since the configurations (|2.4p depend 
only on the dual tori, the original torus seems virtual. It is introduced just for treating all the 
dual tori equally. One could regard one of the dual tori as an original torus. 
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The fermionic matrix ^ is similarly decomposed into blocks as 

(g) V^ii ••• (^i'^(g^i'^\ 



yif^^ (g) tp^^ ■ ■ ■ ip'^^ (g) iIj^'^j 

where ip""^ and ip"'^ represent spinor fields on and T^, respectively. Each 
block ip"''' (g •i/^"^ is in a bifundamental representation {p'^,p^) under the gauge 
group U(p") X U(p*). It turns out [5] that -0"^ has a topological charge on 
as 

3 3 3 / 

1=1 1=1 1=1 ^ 

Indeed, by defining an overlap-Dirac operator, which satisfies a Ginsparg- Wilson 
relation and an index theorenj^l, the Dirac index, i.e., the difference between the 
numbers of chiral zero modes, was shown to take the corresponding value^. In 
the present paper, we do not specify forms of the Dirac operator, and just as- 
sume that the correct number of chiral zero modes arises in the large- limit 
or in the low energy effective theory. 

3 Matrix configurations for some phenomenological 
models without the Higgsino fields 

We now study matrix configurations that provide matter content of some phe- 
nomenological models. 

As we saw in the previous section, in the present formulation, one can realize 
models with a gauge group U(p^) xU(p^) x • • • xU(p'*) and bifundamental matter 
fields. Unfortunately, as we showed in ref. [9], in the model with the gauge group 
U(3) X U(2) X U(l)2 ~ SU(3) X SU(2) x U(l)^ C U(7), either the right-handed 
singlet neutrino or the B-L gauge field cannot be included in the model. Neither 
do we have any solution for the fluxes that provide the SM fermions with three 
generations, in any models whose gauge group is a subgroup of U(7) or smaller 
groups. 

We then consider models whose gauge group is a subgroup of U(8). All the 
subgroups of U(8), which yield the SM fermions with three generations, are 
U(4) X Vl{2) X U^j(2), U,(3) x U(l) x Ui(2) x U^(2), U(4) x Ul(2) x U(1)2, 
and Uc(3) x U(l) x Ul(2) x U(1)^. We will examine all of them. Since the 



(2i 



U^nlmf - mint) = l[{qf - q\) = H (-^{nf - )) . (2.9) 



These techniques were developed in the lattice gauge theories [33] and applied to MM 
and NC geometries [34] , 

The same results were obtained in the fuzzy spheres [351 113| . 
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Uc(3) X U(l) X Ul(2) X U(l)^ model turns out to have no solution, we further 
study a model with Uc(3) x Ul(2) x 11(1)"^ C U(9), in order to find a model 
where the extra gauge group is Abelian. 

As for the Higgs field, the bosonic matrices in the extra dimensions Vi, i.e., 
the gauge fields in the extra dimensions, give scalar fields on our spacetime. 
When the matrices are decomposed into the blocks, some of the block elements 
have the same representation under the gauge group as the Higgs field. However, 
the off-diagonal blocks give massive fields in general. A more serious problem is 
that, as we mentioned in the introduction, even if scalar fields are massless at 
the tree level, it is difficult to keep them massless against quantum corrections, 
which is well-known as the naturalness or the hierarchy problem. In the gauge- 
Higgs unifications [36j, higher-dimensional gauge symmetries protect the scalar 
mass from the quadratic divergences of the cutoff order, but it still can receive 
quantum corrections of the order of the Kaluza-Klein scale (see also ref. |37]). 

In this section, we will find matrix configurations that provide the gauge 
fields and the SM fermions, assuming that the above-mentioned Higgs candi- 
dates or some other mechanisms cause the electroweak symmetry breaking. 



3.1 U(4) X Ul(2) X Uk(2) model 

We first consider the case with the number of blocks being three, i.e., /i = 3, 
and the integers p'^ taken to be 4, 2, 2 for a = 1, ... , h. The gauge group is then 
U(4) X U(2) X U(2). 

The fermionic species are embedded in the fermionic matrix ip as 





Q 


u d 





I 






V e 







o 






o 



where q denotes the quark doublets, / the lepton doublets, u and d the quark 
singlets, and v and e the lepton singlets. The entries denoted as o give no 
massless fermions since, as we will see below, they are set to have a vanishing 
index. The lower triangle part can be obtained from the upper part by the 
charge conjugation transformation. 

The fields q and / are in the (4,2, 1) representation under the gauge group 
U(4) X Ul(2) X Uk(2). The fields n, d, u, and e are in (4,1,2). As we will 
see in subsection 13. 4| the fermionic fields have the correct representation under 
the SM gauge group SUc(3) x SUl(2) x Uy(l), which is a subgroup of U(4) x 
Ul(2) X Uk(2). The color group SUc(3) and the lepton number U(l) are unified 
to SU(4), which is reminiscent of the Pati-Salam model [38j. 
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We now determine the integers qf specifying the magnetic fluxes. From 
(j2.9p . only the differences q^ — q^ are relevant to the topology for the block ip""^ . 
We thus define 



Qi 



„ab 



n 



(3.2) 
(3.3) 



In order for (j3.1|) to have the correct generation number, the integers q""^ must 
have the values 

^0 -3 3^ 

I . (3.4) 

0/ 

The lower triangle part is obtained from the upper part by the relation q°'^ = 
—q^"". The block component with a vanishing index gives no chiral zero modes, 
and thus no massless fermions on our spacetime. Even when there are right- 
handed and left-handed zero modes in a pair, they will obtain a mass through 
quantum corrections, though they are massless at the tree level. 

The fluxes in each T^, g"*, are obtained by solving eq. (j3.3p with the 
fluxes (|3.4p . Here, we note that eq. (|3.3p is invariant under the permutations 
and qf', and also under the sign flips: qf' 



among q'f , 



^ab ^ab 
HI 1 H2 



^ab 
12 ' 



ab 
ab 



^ab. ^ab 

^ab 



„ab 



^ab V ^ab ^ab v ^ab. ^ab , „afo „ab ^ ^ab 

H2 ~^ '/2 > ys ~^ ~% ' yi ~^ HI J 12 ~^ ~H2 ■> 



— )• — . By using these symmetries, we can fix the order of , q2 , and 
'""^ and their overall signs. Under this constraint, there are four solutions for 
""^ We list them in Table [H In order to save space, we have omitted the 



diagonal elements that always vanish. We will call those matrices q^^^ . 
3.2 Uc(3) X U(l) X Ul(2) x U/j(2) model 

We then consider the case where the gauge group U(4) in the previous subsection 
is broken down to Uc(3) x U(l) by the fluxes, i.e., the case with /i = 4 and 
p'^ = (3,1,2,2). 

The fermionic species are embedded in the fermionic matrix as 



o 





q 


u d 




o 


I 


V e 






o 


o 








o 



(3.5) 
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Table 1: Fluxes in each in the U(4) x Ul(2) x \Jr{2) model, i.e., all the 
solutions for eq. (j3.3p with the fluxes ()3.4p . The diagonal elements are 
omitted. 



The fluxes in must have the values 

/o 





\ 



3\ 

3 





(3.6) 



The fluxes in each T^, qf^, are obtained by solving eq. (j3.3p with the fluxes 
(|3.6p . They can be obtained by doubling the first row of the matrices in Table 
dJ However, if qj^ = for all I, which is equivalent to qf = qf for all I, the first 
diagonal block and the second one in the bosonic configuration (j2.4p become 
identical, and the gauge group is enhanced from U(3) x U(f) to U(4), which 
brings us back to the case in the previous subsection. Hence, we must find a 
solution that has both qf'^ = and qf^ ^ 0, depending on I. There are four 
solutions, which we list in Table [2j 



3.3 U(4) X Ul(2) X U(l)2 model 

We next consider the case where the gauge group Ur(2) is broken down to 
U(l) X U(l) by the fluxes, i.e., the case with h = A and p"- = (4, 2, 1, 1). 
The fermionic species are embedded in the fermionic matrix ip as 

( 





Q 


u 


d 





I 








u 


e 

































(3.7) 
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/o 
V 



/o 



V 



qf 



-l\ 
1 

0/ 



qf 



(-2 



(2 1 



/4 



(2 



\ 




c-ab 
I3 



/o 




/6 3 







/4 1 



-3 
0/ 



/-2 



V 




Table 2: Fluxes in each in the Uc(3) x U(l) x Ul(2) x \Jr{2) model, i.e., ah 
the solutions for eq. (j3.3p with the fluxes (13. 6p . The diagonal elements are 
omitted. 



The fluxes in must have the values 

/O -3 







V 



3 3\ 





(3.8) 



The fluxes in each T^, qf, are obtained by solving eq. (|3.3|) with the fluxes 
(|3.8p . They can be obtained by suitably doubling the last column of the matrices 
in Table [TJ Under the constraints mentioned in the previous subsections, there 
are nine solutions, where we have also fixed the exchange between the last two 
columns. We list them in Table O 



3.4 Uc(3) X U(l) X Ul(2) x V{lf model 

We now consider the case where both U(4) Uc(3) x U(l) and U_r(2) U(l)^ 
take place by the fluxes, i.e., the case with h = 5 and = (3, 1, 2, 1, 1). 
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qf 


qf 


qf 




/-I 1 -l\ 

2 
V -2/ 




{ 


'l 1 -l' 
-2 


) 




/3 3 3\ 


V 0/ 






/-I 1 -l\ 

2 
V -2/ 






/l 1 l\ 


V 0/ 






/3 3 -3\ 
-6 

V -6/ 






r = i) 






/I 1 l\ 



V 0/ 






/3 3 1 \ 
-2 

V -2/ 






/-I 1 -3\ 
2 -2 

V -4/ 






/l 1 l\ 



V 0/ 






/3 3 -l\ 

-4 

V -4/ 






t • 3 






/I 1 l\ 



V 0/ 




( 


3-3 1 
-6 -2 
4, 


) 


1 
\ 


-1 -3 
-2 
-2 


) 




/i 1 i\ 



V 0/ 




/ 
\ 


'3 -3 -1^ 
-6 -4 
2, 


/ 




( • 3 




\ 


1-3 1 
-4 


) 


/ 
I 


3 1 1 
-2 -2 
0, 


) 


\ 


-1 -1 




) 


( 


'l -3 3' 
-4 2 


) 


( 


'3 1 -1^ 
-2 -4 

-2, 


) 


( 


'-1 -1 -3 
-2 


) 


/ 
\ 


'13 1^ 
2 

-2. 


) 


( 


'3 -1 -1^ 
-4 -4 
0, 


) 



Table 3: Fluxes in each in the U(4) x Ul(2) x U(1) x U(1) model, i.e., ah 
the solutions for eq. ()3.3p with the fluxes (13. 8p . The diagonal elements are 
omitted. 
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The fermionic species are embedded in the fermionic matrix ^ as 



/ o 



V 



(3.9) 



o J 



The fluxes in must have the values 



/o 


-3 3 


3\ 





-3 3 


3 


















v 




0/ 



(3.10) 



The fluxes in each T^, qf^ , are obtained by solving eq. (|3.3p with the fluxes 
(|3.10p . They could be obtained by suitably doubling the last column of the 
matrices in Table [21 or by suitably doubling the first row of the matrices in 
Table [3l Unfortunately, there is no solution, however. 
We then generalize the fermion embedding (|3.9p to 



(o 



u 
o 



u 
v(v) 
V 
o 



\ 

(o 



o 
e' 



(3.11) 



\ 



\ 



(3.12) 



u q d\ d2 

O I ei 62 

o o o 
o i'{r') 

o I 

though they cannot be applied to the models with the larger gauge groups in 
the previous subsections. Here, we have omitted the lines that separate the 
matrix elements since they are no longer necessary. 

In fact, (j3.1ip and (|3.12p are the most general embeddings where all the 
elements have the correct representation under the SM gauge group SUc(3) x 
SU2,(2) X Uy(l): One can easily see that they have the correct representation 
under the SUc(3) x SUl(2). Since we have five U(l) gauge groups coming from 
each diagonal block, we consider their linear combinations 

5 

(3.13) 

i=l 



12 





Q 


u 


u' 


d 


I 


/' 




e 


e' 


Y 


1/6 


2/3 


2/3 


-1/3 


-1/2 


-1/2 





-1 


-1 


B 


1/3 


1/3 


1/3 


1/3 

















L' 








-1 





1 





1 


1 





Ql 


1 











1 


1 











Q'r 





1 





1 







-1 




1 


1 





















x^ 










Y 


1/6 


-1/2 







1/2 


1/2 










B 


1/3 
























L' 





1 





















Ql 








-1 


















Q'r 











-1 




-1 







Table 4: U(l) charges for for g, n, n', d, /, vif)-, e, and e' in ()3.1ip . and the 
corresponding values for the coefficients in ()3.13p . 

where Q^ is the U(l) charge from the ith block. For the (j3.1ip case, we can 
consider the hypercharge Y , baryon number B, lepton number L', left-handed 
charge Ql, and right-handed charge Q'pi^. Their charge for q, n, u', d, /, /', viv), 
e, and e' , and the corresponding values for are given in Table [H The charges 
L' and Q'j^ have the proper interpretation when u' , I' , and e' disappear, and ly, 
not v', is chosen in Tabled and thus in (|3.1ip . This case reduces to the original 
embedding ()3.9p . This also ensures that the fermion species in (j3.ip . ()3.5p . and 
p.7p in the previous subsections have the correct representation under the SM 
gauge group and the extra U(l) gauge groups, which are subgroups of the gauge 
group in the corresponding model. 

For the (j3.12p case, we consider the hypercharge Y, baryon number B, left- 
handed charge Ql, and other two charges Q'l and Q2- Their charge for q, u, di, 
d2, I, ei, and 62, and the corresponding values for are given in TableO 
The charges Q'^ and Q'2 have no proper interpretations. We also note that both 
in the cases (j3.1ip and ()3.12p . a linear combination of those five U(l) charges 
gives an overall U(l) and does not couple to the matter. Only four U(l) charges 
couple to the matter. 
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Q 


u 


di 


d2 


I 




ei 


62 


Y 


1/6 


2/3 


-1/3 


-1/3 


-1/2 





-1 


-1 


B 


1/3 


1/3 


1/3 


1/3 














Ql 


1 











1 











Q'l 








1 











1 





Q'l 











1 





1 





1 



Y 



B 



Ql 



Q'l 



Q'2 



1/6 



1/3 







-1/2 











1/2 







1/2 



Table 5: U(l) charges for for u, di, d2, I, vi^)-, ei, and 62, in (I3.12p . and the 
corresponding values for the coefficients in (|3.13p . 



The fluxes in T must have the values 



/O X 


—3 3 — X 


3 


\ 





y-3 ±3 


3 - 


z 




y 












z 




V 







J 



for (IXTTD . and 



/O 



-3 
-3 



3 
3 



x\ 

y 





±3 




(3.14) 



(3.15) 



for (j3.12p . Here, x, y, and z can take an integer 0, 1, 2, or 3. Since the up-type 
quarks are embedded at the two places u and u' in (j3.1ip . the corresponding 
fluxes can take several values x and 3 — x in (I3.14|) . The same is true for 
/ and and so on. The double sings of ±3 are chosen whether v or is 
embedded in (j3.1ip and (|3.12p . The fluxes in each T^, q^'', are obtained by 
solving eq. ([33]) with the fluxes (|3Tl]) and ([313]) . As we will show in 
Appendix lA. 21 there is no solution, either. We thus conclude that there is no 
solution in the Uc(3) x U(l) x Ul(2) x U(1)2 model. 
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3.5 Uc(3) X Ul(2) X U(l)^ model 



We then consider the case with the gauge group Uc(3) x 11^(2) x U(l)^ C U(9), 
i.e., the case with h = 6 and p'^ = (3, 2, 1, 1, 1, 1), in order to find a model with 
the SM gauge group plus extra U(l)'s. 

The fermionic species are embedded in the fermionic matrix ip as 
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(3.16) 



(3.17) 



(3.18) 



where v can be either v or v. They exhaust all the embeddings that have the 
correct representation under the SM gauge group. 

The fluxes in T^, q"^ can take several values for (j3.16p . (|3.17|) . and (j3.18p . 
as in (j3.14p and (j3.15p . We can also determine the fluxes for T^, g^'', by solving 
eq. (j3.3p . As we will show in Appendix lA. 31 there are fifteen solutions for the 
case ([3l6]) . The fluxes are listed in Table [H Table HZl and Table [13 The 
corresponding fluxes are written in (jA.lSp to (|A.20p . 

For the case (j3.17p . there are two solutions. The fluxes are listed in 
Table [T5j The corresponding fluxes are (|A.10p and (jA.lip . For the case 
p.lSp . there is no solution. 
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4 Matrix configurations for the phenomenological mod- 
els with the Higgsino fields 

In this section, we assume situations where the Higgs mass is protected by the 
supersymmetry possessed by the IIB MM somehow. Then, we wih find matrix 
configurations that yield candidates for the Higgsino fields as well as the gauge 
fields and the SM fermions. We just find candidates for fermion species, and 
do not study how the surpersymmetry is realized in the whole spectrum. 



4.1 U(4) X Ul(2) X Ur{2) model 

We first consider the model studied in subsection 13.11 The Higgsino fields h^, 
and hd can be added to the fermion embedding (j3.ip 
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u d 
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V e 







hu hd 






o 



The fluxes must be altered from ()3.4p to 

^0 




(4.1) 



(4.2) 



with ^ 0. There are three solutions for the fluxes. They are listed 
in Table [6l where we also write the values of . Unfortunately, there is no 
solution with q^ = 1. 



4.2 Uc(3) X U(l) X Ul(2) x Ur(2) model 

We then consider the model studied in subsection 
and hd are added to the embedding (|3.5p as 
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hu hd 
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The Higgsino fields h^ 



(4.3) 



A similar embedding was given in ref. 
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Table 6: Fluxes in each in the U(4) x 11^(2) x \Jii{2) model with the Higgsino 
candidates, i.e., all the solutions for the fluxes (j4.2p . The diagonal elements 
are omitted. The flux for the Higgsinos, , is also listed. 



The fluxes in have the values 



/o 


-3 3 \ 





-3 3 







V 


0/ 



(4.4) 



with q^ 7^ 0. Solutions for the fluxes can be obtained by suitably doubling 
the first row of the matrices in Table [H There are four solutions. They are 
listed in Table [71 where we also write the values of q^ . 

Because the lepton doublet / and the Higgsino hu in (j4.3p have the same 
representation under the SM gauge group, they can substitute for each other. 
Then, the fluxes can have more general values 



/o 



v 



-3 



3\ 



3 




(4.5) 



Solutions for the fluxes can also be obtained by using the the matrices in 
Table [6l There are fourteen solutions, other than those in Table [71 We do not 
present them here, in order to save space. 
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32 



Table 7: Fluxes in each in the Uc(3) x U(l) x Ul(2) x Uh(2) model with the 
Higgsino candidates, i.e., all the solutions for the fluxes (j4.4p . The diagonal 
elements are omitted. The flux for the Higgsinos, , is also listed. 



4.3 U(4) X Ul(2) X U(l)2 model 

We next consider the model studied in subsection 13.31 The Higgsino fields h. 
and hd are added to the embedding (13. 7p as 





Q 


u 


d \ 





I 








V 


e 







hu 




















) 



(4.6) 



The fluxes in must have the values 



/O -3 
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3 \ 
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/ 



(4.7) 



with q^^,q^'' ^ 0. Solutions for the fluxes can be obtained by suitably 
doubling the last column of the matrices in Table [6l There are five solutions. 
They are listed in Table [HI where we also write the values of q^"^ and g^'*. 
Solutions with q^^ = and q^'^ ^ 0, or q^^ ^ and q^''^ = 0, can also be 
obtained by suitably combining the last column of the matrices in Table [T] and 
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Table 8: Fluxes in each in the U(4) x Ul(2) x U(1)^ model with the Higgsino 
candidates, i.e., solutions for the fluxes (14. 7p (Part 1). The diagonal elements 
are omitted. Also listed are the fluxes for the Higgsinos, g^" and q^'^, both 
of which take nonzero values. 



Table [6l There are eight solutions. They are listed in Table O where we also 
write the values of g^" and g^'*. 

4.4 Uc(3) X U(l) X Ul(2) x \]{lf model 

We now consider the model studied in subsection 13.41 The Higgsino fields hu 
and hd are added to the embedding (|3.9p as 
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(4.8) 
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Table 9: Fluxes in each in the U(4) x Ul(2) x U(1)^ model with the Higgsino 
candidates, i.e., solutions for the fluxes (14. 7p (Part 2). The diagonal elements 
are omitted. Also listed are the fluxes for the Higgsinos, q^'^ and q^''-, one 
of which takes nonzero values. 
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The fluxes in have the values 



/o 


-3 3 


3 \ 





-3 3 
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(4.9) 



with q^^,q^'^ ^ 0. Solutions for the fluxes can be obtained by suitably 
doubling the last column of the matrices in Table [71 or by suitably doubhng 
the first row of the matrices in Table [HI There is one solution. It is listed in the 
first row of Table [TOl where we also write the values of q^"^ and q^'^ . Solutions 
with q^" = and q^'^ ^ 0, or g^" ^ and q^'^ = 0, can also be obtained by 
suitably combining the last column of the matrices in Table [2] and Table [Tj or 
by suitably doubling the first row of the matrices in Table [9j There are two 
solutions, which are listed in the second and the third rows of Table [TOl 

As in the Uc(3) x U(l) x 11^^(2) x U/j(2) model, since I and hu have the same 
representation under the SM gauge group, the fluxes can have more general 
values 

/O -3 



q^ 
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±3 
qHu 
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3 \ 
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/ 



(4.10) 



Solutions for the fluxes can be obtained by using the matrices in Table [8] 
and Table [H There are three solutions. We list them in Table [Til where we 
also write the values of q^"^, q^'^ and q^. As for the double sign of ±3 in (j4.1Up . 
all of the solutions take +3. 



5 Probability distribution over the phenomenologi- 
cal models 

We now study the dynamics of MM semiclassically, and estimate the probabili- 
ties for the appearance of the phenomenological models obtained in the previous 
sections. 



5.1 Model and instanton actions 

We first specify the action of the MM. We here consider a ten-dimensional 
NC torus with an anisotropy of sizes between four and six dimensions, while 
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Table 10: Fluxes in each in the Uc(3) x U(l) x Ul(2) x \J{lf model with 
the Higgsino candidates. All the solutions for the fluxes (j4.9p . The diagonal 
elements are omitted. The fluxes for the Higgsinos, q^" and g^", are also 
listed. 
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Table 11: Fluxes in each in the Uc(3) x U(l) x Ul(2) x \J{lf model with 
the Higgsino candidates. All the solutions for the fluxes (j4.10p . other than 
those listed in Table [TOl The diagonal elements are omitted. The fluxes for 
the Higgsinos and the lepton doublet, g^", and q^, are also listed. 
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one can alternatively consider a six-dimensional torus with a four-dimensional 
uncompactified spacetime. The bosonic action is given by the twisted Eguchi- 
Kawai model |26t I27j (see, for instance, ref. |28J) as 

Sb = -!3M ^oi tr (v. V, V\ Vj) - /3'A/- ^ Z,^ tr (v^ V. V^) 

-/3"M J2 tr (v. V] V^) + Z,^ tr (v^ V, v/)] , (5.1) 

where and Vj are U(A/') matrices with //, = 0, . . . , 3 and i, j = 4, . . . , 9. We 
take twists 

Zoi = exp (27ri^) , ^23 = exp (27ri^) , 

iV4 iVs 
2:45 = exp (27ri-^) , Zqj = exp (2vri-^) , Zgg = exp (27ri-^) ,(5.2) 

with a huge difference in Ni between / = 4, 5 and / = 1, 2, 3: 

A^4 ^ iVs > A^i - =^ iVs • (5.3) 

The other twists are taken to be one. The total matrix size M is taken to be 

5 

M = k\{Ni. (5.4) 

1=1 

Note that the matrix size (|5.4p is k times larger than is usually expected from 
the integers that specify the twists ()5.2p . 

We now consider background configurations as 

V^ = \(^Vi, (5.5) 

where are U(A'4A'^5) matrices and Vi are \]{kNiN2N^) matrices. The size of 
our spacetime and that of the extra dimensions are given by eA'^; with / = 4, 5 
and I = 1,2, 3, respectively, where e is a lattice spacing. They thus have a huge 
anisotropy due to (|5.3p . 

We then consider the topological configurations ()2.4p for Vi. In fact, they 
are classical solutions for the action ()5.ip (see, for instance, ref. |40)). In order 
to match the matrix size, 

h 3 
Yn>-2ntp'' = kl{Ni (5.6) 

a=l 1=1 

is required. Plugging ()2.4p into ()5.ip . we obtain the classical action 

^ ' ' si mf 

1=1 a=l 



Sb = -2f3J^N^N, Y E <"2^3P" cos i^27T^ji- + -L)) , (5.7) 
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where we have written only the first term in (jS.ip . 

Recall that the integers and are specified by an original torus with 
the integers Ni and si via (|2.6|) or (|2.7p . We now consider the configurations 
whose integers A'^; and s/ agree with the ones that specify the twists Zij in the 
action (jS.ip via ()5.2p . The other configurations have much larger actions. Then, 
from (j2.6p and ()2.7p . we can find the relation 

si mf qf ^ f ^ ^ 

N'l^'nf ~ Ninf ~ ~2r [Wi ~ n° 

By plugging (15. Sp into (j5.7p . we find that the classical action (15. 7p takes the 
minimum value if and only if 

q1 = 0^nf = Ni (5.9) 
for Va and V/. Then, the constraint (15. 6p becomes 

h 

^p'' = k. (5.10) 

a=l 

Therefore, if we choose the parameters of the MM action (j5.ip as in (j5.4p . block 
diagonal configurations, where the total number of the blocks is specified by 
(jS.lOp . are dynamically favored. 

We then consider small fluctuations around the minimum, i.e., configura- 
tions with \qf\ <^ Ni. The classical action (15. 7p is approximated as 

AS, 4^'(3^ f: 7^. iZp'^ilf)' > (5-11) 

1=1 ^ '■' a=l 

where we have written the difference from the minimum value. We hereafter 
will call it an instanton action since it is a classical action of a topological 
configuration. 

For comparison, let us consider cases with large fluctuations, for instance, 
configurations where the total number of blocks is different from (j5.10p . and 
configurations specified by original tori with the integers A'^;, si which are dif- 
ferent from the twists in the action (j5.ip . In those cases, the classical action 
(j5.7p receives an enhancement factor of order (Ni)'^, compared to (j5.1ip . 

5.2 How to take large- limits and the probability distribution 
over the string vacuum space 

We now consider relations between the prescription of how to take a large-N 
limit in the MM and the probability distribution over various matrix configura- 
tions, i.e., various string vacua, based on the semiclassical analysis. For details, 
see ref. [9j. 
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We first compare the IIB MM action (j2.ip and the unitary MM action 
(j5.ip . by considering a correspondence between the Hermitian matrices and the 
unitary matrices as 

Vm ~ exp ( 27ri^ ) . (5.12) 



eNi 

By plugging (j5.12p into (jS.ip . and comparing it with ()2.ip . we find a relation 
among the coupling constants in (j5.ip and (j2.ip as 

2 VNiJ SiiBMM 

with / = 1,2,3. Similar relations can be obtained for /3' and /3". 

Since both ()5.1ip and ()5.13p depend on f3/{Ni)^, by defining a combination 

5lIBMM _ 1 



the instanton action (|5.1ip becomes 

3 h 



(5.14) 



«=1 a=l 



It then follows that scaling limits of fixing Oiibmm-^"/^^ with a>— l,a = — 1, 
and a < —1 make the coefficient of the instanton action (|5.15p infinite, finite, 
and vanishing values, respectively. Together with fixing the torus size eM^^^, 
those scaling limits correspond to fixing fl'nBMM-^'^ with 7 = a + 4/5. The three 
cases give drastically different results: 

1. If we take a large- limit by fixing ffnBMM-^"/^^ with a > —1, or by 
fixing ^iiBMM-^''^ with 7 > —1/5, the instanton action (j5.15p diverges for 
nonzero qf, and thus only a single topological sector survives. While in 
the present model, the topologically trivial sector, qf = 0, is chosen, in 
more elaborated models, desirable sectors may be chosen uniquely by the 
dynamics. This is drastically different from the situations where physicists 
usually consider the landscape. 

2. In a limit with a<— lor7<— 1/5, the instanton action (|5.15p vanishes 
for all qf, and all the topological sectors appear with equal probabili- 
ties. Then, the estimation for the probability distribution over the string 
vacuum space reduces to the number counting of the classical solutions. 
Moreover, in a limit with a < —1 — 2/5, a still larger number of config- 
urations, which were studied as the large fluctuations below ()5.1ip . can 
also appear. 
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3. In a limit with a = —1 or 7 = —1/5, the instanton action (j5.15p takes 
finite values for finite qf, and all the topological sectors appear with finite 
and different probabilities. 

5.3 Probabilities for the appearance of the phenomenological 
models 

We now estimate the probabilities for the appearance of the phenomenological 
models obtained in sections [3] and HI 

We first consider the U(4) x U_l(2) x Ur(2) model with the fluxes given by 
the first row in Table [TJ By solving (j3.2p . qf are determined as 

Ql = {qi,qi + l,qi-l) , 

92 = (92,92 - 1,92 - 1) , 

93 = (93,93 - 3,93-3), (5.16) 

for a = 1, . . . , /i = 3. Since only the differences are specified in (j3.2p . qf are de- 
termined with arbitrary integer shifts gi, 92, and 93. The instanton action (j5.15p 
takes various values depending on those arbitrary integers qi. The minimum 
value 

AS, = ^28 , (5.17) 

is attained by gi = 0, 92 = or 1, and 93 = 1 or 2. 

One could further lower the instanton action (|5.15p to 

AS, = ^24 , (5^18) 

by choosing the fractional values qi = 0, 92 = 1/2, and 93 = 3/2. This corre- 
sponds to modifying the twists in the action (j5.ip from (j5.2p to 

^45 = exp ^27ri-^ 



N2 2iV| 

S3 1 



exp ( 27r^ ( ^ + — ^ ) ) , (5.19) 



although they are not natural in the present MM. 

The probability of the appearance of this model is semiclassically given as 
g-A56^ multiplied by a factor coming from quantum corrections. There exist 
configurations with the same action, but with and qf different from (j5.16p . 
and thus the probability must be divided by this numerical factor. 
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In the same way, we estimate the minimum instanton action for the other 
cases in Tabled) We also examine the Uc(3) x U(l) x 11^(2) x IJr(2) model with 
the fluxes given in Tabled and the U(4) x Ul(2) x U(1)^ model with the fluxes 
in Table [31 The results are drawn in Table [T2l We further examine the models 
with the Higgsino fields, studied in section [H and list the results in Table [13l 
We do not present the results from the fluxes (|4.5p and (|4.10p . since they 
give larger values of the instanton actions. 

Within these models and in the integral qi case, the minimum instanton 
action ^^^^-ASj, takes various values between 28 and 44 over the various models. 
The minimum of the minimum values, 28, is attained by the first solution in the 
U(4) X Ul(2) X Ur(2) model, the last one in the Uc(3) x U(l) x Ul(2) x Vr{2) 
model, and the first one in the U(4) x U2,(2) x U(l)^ model in Table [T2| and also 
by the second solution in the U(4) x U/, (2) x U/j(2) model and the seventh one in 
the U(4) X Ul(2) x U(1)^ model in Table [T3l We cannot find drastic differences 
among the models with various gauge groups, and between the models with 
and without the Higgsinos. 

There also exist solutions with lower instanton actions, which do not yield 
the SM fermions. Moreover, as far as we consider the models whose gauge 
group is a subgroup of U(8), i.e., the MM (15. ip with the parameter k = 8, the 
U(8) model without any gauge symmetry breaking has the lowest action. In 
order to make phenomenologically-attractive models most probable, we need 
to elaborately modify the MM (j5.ip . However, we can also find that the phe- 
nomenological models have rather small instanton actions and are sufficiently 
probable if the coefficient A is not so large. 

6 Conclusions and discussions 

In this paper, we considered the situations where the IIB MM is compactified 
on a torus with magnetic fluxes, and exhausted all the matrix configurations 
that yield all the phenomenological models whose gauge group is a subgroup 
of U(8), with and without the Higgsino fields. We then studied the dynam- 
ics of MM semiclassically, and estimated the probability distribution over the 
phenomenological models. 

There remain some important problems. While we found the embedding of 
the Higgs field in the matrices, we need to consider how the electroweak sym- 
metry breaking occurs. The gauge-Higgs unifications and the recombinations 
of the intersecting D-branes (see, for instance, ref. |41) ) are close to the present 
situation and may be helpful to consider this issue. We should also study val- 
ues of the Yukawa couplings and the flavor structure. They can be obtained by 
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gauge group 


integral qi 


fractional qi 


U(8) 








U(4) X Ul(2) X VRi2) 


28 


24 




44 


40 




36 


32 




36 


32 


Ue(3) X U(l) X Ul(2) X \Jr{2) 


44 


43 




40 


39 




36 


36 




28 


27 


U(4) X Ul(2) X U(l)2 


28 


25 




40 


37 




32 


29 




36 


32 




44 


40 




40 


37 




36 


36 




40 


37 




36 


33 



Table 12: Values of the instanton actions ^^-^ • ASb in the U(8) model, the 
U(4) X Ul(2) X Vr{2) model, the Uc(3) x U(l) x Ul(2) x U/j(2) model, and the 
U(4) X Ul(2) X U(l)^ model, without the Higgsino fields. Each row of this table 
corresponds to those in Table [H Table [21 and Table [3l The second column 
shows the minimum value within the integer values of qi as in (j5.17p . while 
the third column is the results within the fractional values of qi as in (j5.18p . 
Fractional values of qi are not natural in the present MM. 
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0, -8 
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0, 32 


40 


37 
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Ue(3) X U(l) X Ul(2) X U(l)2 


-8, -8 


36 
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0, -8 


40 


39 




0, 32 


40 


39 



Table 13: Values of the instanton actions • ASb in the the U(4) x Ul(2) x 
Ur(2) model, the Uc(3) x U(l) x Ul(2) x Vr{2) model, the U(4) x Ul(2) x 
U(l)2 model, and the Uc(3) x U(l) x Ul(2) x 11(1)^ model, with the Higgsino 
candidates. Each row of this table corresponds to those in Table El Table O 
Table [HI Tabled and Table [TOl We do not present the results from the fluxes 
(j4.5p and (|4.10p . since they give larger values of the instanton actions. The 
fluxes for the Higgsinos, q^" and q^"^, are also listed. The third column shows 
the minimum value within the integer values of qi as in (|5.17p , while the fourth 
column is the results within the fractional values of qi as in (I5.18p . Fractional 
values of qi are not natural in the present MM. 
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calculating overlaps among the zero-mode fields. Related works were given, for 
instance, in ref. [i2t |39]. 

While we considered situations where the supersymmetry protects the Higgs 
mass, we need to study how to keep a part of the supersymmetry possessed by 
the IIB MM and how to break it at low energies. It is also possible that 
the supersymmetry is broken at high energies but leave some remnants (see, 
for instance, ref. [12]). Other resolutions of the naturalness or the hierarchy 
problem are also interesting to consider, such as composite models, models 
with large-extra dimensions, some stringy or quantum-gravitational effects. 

The models we considered in the present paper have extra U(l) gauge groups 
and are anomalous within the gauge dynamics. This anomaly may be canceled 
via the Green-Schwarz mechanism by the exchange of RR-fields, which can also 
make the extra gauge fields massive. In order to realize this, some constraints 
as the RR tadpole cancelation condition should be imposed on MM. We should 
also generalize topological configurations by introducing the complex structures 
and the Wilson lines. 

While we assumed toroidal compactifications and worked in a unitary matrix 
formulation in this paper, we should study relation between the unitary MM and 
the IIB MM, and how the parameters in the unitary MM action are determined 
from the IIB MM. We also need to consider how to interpret spacetime and 
matter in the matrices and how to describe compactifications in the matrices. 

We will come back to these issues in future publications. We hope that 
these studies will give us some helps for both exploring the phenomenological 
models and studying the formulations of MM. 
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A Solutions of the fluxes qf 

In this appendix, we find solutions of the fluxes qf^ by solving eq. (|3.3p in 
the Uc(3) X U(l) X Ul(2) x \]{lf model and the Uc(3) x Ul(2) x U(1)^ model 
without the Higgsino fields, studied in section [3l 
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A.l Two theorems 



Before we solve eq. (j3.3p in the concrete models, we prove two theorems. 

Theorem 1: Pick up three elements from the fluxes and focus on g"^, 
q^'^, and q"''^: 

:::) ■ 

Then, it is impossible that all of Iq""^], {q^'^l, and \q"''^\ take 1, 2, or 3. In other 
words, if |^"^|, [q^'^l, and \q"''^\ are 0, 1, 2, or 3, then, at least one of them must 
be 0. There is one exception: [g^'^l, \q°''^\) = (2,2,2) is possible. 

Proof: Let us first consider the case |(?^^|) = (1,1)- Then, the 

fluxes \qf''\ must take 1 for all I. So must Igf^^l. It then follows from qf' 



Ml 



that must be or 2. Then, \q'^'^\ is or 8. Hence, 
(1, 1, 1), (1, 1, 2), or (1, 1, 3) is not possible. 

Similarly, by considering the cases \q^'^\) = (1,2), (1,3), (2,2) 

and (3,3), one can prove the theorem. 

Theorem 2: Consider the fluxes 



1) = 
(2,3), 



which is equivalent to 



/o 



V 

/O - 



V 



3\ 

X 





3\ 
3 



(A.2) 







X 



(A.3) 



as can be seen by exchanging the second and third rows and columns. Then, 
it is impossible that x takes 0, 1, ±2, or —3. Within |x| < 3, only x = —1 and 
X = 3 are possible. 



Proof: For 




(A.4) 



the solutions for eq. 
for 



are listed in Table [T] One can also list the solutions 




0. 



(A.5) 
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qf 



/-2 -1 



/O 



V 



Q2 



/2 



V 




/-2 1 
3 

V 



/-2 1 
3 

V 



/O 



-aft 
93 

3 
3 



/2 




Table 14: Fluxes in each for the fluxes fO]) with x 
elements are omitted. 



-1. The diagonal 



By combining two solutions from each, one can construct a solution for ()A.2p . 
There is no solution for the cases x = 0,1, ±2,— 3. For x = 3, the solutions 
are listed in Table [2l For x = —1, there are three solutions, which we list in 
Table m 



A. 2 Fluxes in the Uc(3) x U(l) x Ul(2) x model 

We now study the Uc(3) x U(l) x Ul(2) x U(1)^ model considered in subsection 
13.41 One can immediately find that (I3.15P has no solution for (13. 3p , by applying 
the Theorem 1 to q^'^ = 3, q"^^ = -3, and q^^ = -3. 

For p.l4p . by applying the Theorem 1 to q^^ = 3, one can find that x = 
2; = 0orx = z = 3is allowed. Also, by applying the Theorem 1 to q^^ = —3, 
x = y = 0oTx = y = 3is allowed. Combining these two results, x = y = z = 
or x = y = z = 3 is allowed. They correspond to 
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0/ 
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3\ 










T3 












3 
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V 








0/ 



(A.6) 



(A.7) 
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respectively. In fact, ()A.6P and ()A.7p are equivalent, as one can see by exchang- 
ing the second and the fourth rows and columns, i.e., by exchanging the first 
and the second U(l). 

By applying the Theorem 2 to the first 4x4 part of (IA.6|) , one can find that 
the plus sign in the double sign must be chosen. After all, the problem comes 
back to the case in (j3.10p . Then, we have to conclude again that there is no 
solution in the Uc(3) x U(l) x Ul(2) x U{lf model, even though the fermion 
embedding is generalized to (j3.1ip and ()3.12|) . 



A.3 Fluxes in the Uc(3) x Ul(2) x U(1)^ model 

We then study the Uc(3) x Ul(2) x U(1)^ model considered in subsection 13.51 
One can immediately find that the case (|3.18|) has no solution for q'j^^, by ap- 
plying the Theorem 1 to q^^ = —3, q^^ = 3, and q^"^ = 3. 

For the case (j3.17p . by using the Theorem 1, the possible fluxes turn out 
to be as follows: 



/o 



„ab 
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3 - 






































/o -3 
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o\ 
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\ 



(A.9) 



where x, y, and z can take an integer 0, 1, 2, or 3, while u and v can take 0, 
±1, ±2, or ±3. In (|A.9p . x + y + z = 3 and \u\ + \v\ = 3 must be satisfied. 

For (jA.Sp . by applying the Theorem 2 to the first 4x4 part, i.e., q""^ with 
1 < a, 6 < 4, one can see that the plus sign in the double sign must be chosen. 
Then, g''^ with 1 < a, 6 < 4 are the same as (|3.6p , as one can see by exchanging 
the second and the third rows and columns in g"^. Thus, we can construct 
solutions for (jA.Sp by using the solutions in Table [2j There is one solution. The 

fiuxes are written in the first row in Table [T5j The corresponding fiuxes 
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are (IXSl) with 



1, i.e., 



/o 




3 
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1 
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2 























0/ 



(A.IO) 



v 

For ()A.9p . by applying the Theorem 2 to the first 4x4 part, i.e., q"'^ with 
1 < a, 6 < 4, one can find that u must take —1 or 3. In the same way, x must 
be —1 or 3. As we showed below ()3.10p . there is no solution for {u,x) = (3,3). 
Since x must be positive, the remaining possibility is {u,x) = (—1,3). It then 
follows from x + y + z = 3 and \u\ + \v\ = 3 that y = z = and |?;| = 2. 
One can construct solutions by using those in Table [2] and Table [TH There is 
one solution. The fluxes are written in the second row in Table [151 The 
corresponding fluxes are 
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0/ 



(A.ll) 



For the case (|3.16p . by using the Theorem 1, the possible fluxes turn out 
to be as follows: 
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(A.13) 



34 





~a6 


-aft 
^3 


/-I -2 -1 -1 -l\ 
-10 
1 1 1 


V 0/ 


/-I -2 -1 -1 -l\ 
-10 
1 1 1 


V 0/ 




^-3 3 1 2 ^ 
3 6 4 5 
3 1 2 
-2 -1 
1 






/-I -1 1 -l\ 
10 2 

-1 1 -1 

2 

I -2/ 


/l -2 -1 1 \ 
-3 -2 -1 
13 2 
2 1 

V -1/ 




/3 2 3 3 2 \ 
-10 0-1 
1 1 
-1 

V -V 





Table 15: Fluxes in each in the Uc(3) xUl(2) xU(1)^ model with the fermion 
embedding (j3.17p . The first and the second rows of this table correspond to the 
fluxes (lA.lOp and (lA.lip . respectively. The diagonal elements are omitted. 
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3 — X 




(A.14) 



where x can take an integer 0, 1, 2, or 3, while u, v, and w can take 0, ±1, ±2, or 
±3. In (lAl2ll and (fAl3l) . \u\ + \v\ = 3 is required. In (lAli]) . 
is required. 

Since ()A.12p has a similar structure to (jA.Sp . as can be seen by permuting 
the last three columns, one can find solutions in the same way. There are three 
solutions: The fluxes of the first row in Table [16] yield the fluxes (|A.15p . 
The second and the third rows in Table [16] give ()A.16P . 
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/o 




3 




1 
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(A.16) 



v 







For (|A.13p . q"'^ with a,b = 1, 2, 5, 6 are the same as (j3.8p . and thus one can 
use the results in Table [3l We just simply find solutions for the first 4x4 part 
in (jA.lSp . We then combine the results, and obtain solutions for ()A.13p . There 
are three solutions: The fluxes of the fourth and the fifth rows in Table [16] 
yield fluxes (|AT7l) . The sixth row in Table [16] gives (jA.lSP . 
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1 2 
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3 3\ 
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(A.18) 



For ()A.14p . q with a,b = 1, 2, 3, 4, 6 are the same as the first 5x5 part in 
()A.9p . One can thus follow the same arguments there. There are nine solutions: 
The fluxes in Table [IT] yield the fluxes (IXl9]l . and those in Table [IS] give 
(ICT]) . 

/o -3 3 
3 
-1 
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(A.19) 
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/o 







±2 




V 



3\ 


3 



V 



(A.20) 



In all of the above solutions ([AlO]) . ([AH]) . fO5l) - (|A:20]) . some of the 
fermion species are placed in several matrix elements, i.e., not all of the fermion 
species can be put in a single place with three generations. This is inevitable 
since there is no solution for (j3.10p . It may give some interesting results in the 
flavor structure. 
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Table 16: Fluxes in each in the Uc(3) x Ul(2) x U(1)^ model with the 
fermion embedding (j3.16p (Part 1). The first row of this table corresponds to 
the fluxes (jA.lSp . the second and the third rows to ()A.16p . the fourth and 
the fifth rows to ()A.17p . and the sixth row to (jA.lSp . The diagonal elements 
are omitted. 
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Table 17: Fluxes in each in the Uc(3) x Ul(2) x U(l)'^ model with the fermion 
embedding (j3.16p (Part 2). The solutions in this table correspond to the 
fluxes (jA.lOp . The diagonal elements are omitted. 
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Table 18: Fluxes in each in the Uc(3) xUl(2) xU(l)^ model with the fermion 
embedding (j3.16p (Part 3). The solutions in this table correspond to the 
fluxes ()A.20|) . The diagonal elements are omitted. 
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